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Abstract 

We examine a scenario where the new physics at the LHC includes an approximate 
conformal field theory, where some of the degrees of freedom (aka "unparticles" ) carry a 
color charge. We present a simple argument showing that the production cross section 
for scalar unparticles mediated by a gauge interaction is given by 2 — d times the 
standard particle expression, where d is the scaling dimension of the unparticle field. 
We explicitly check that this is indeed the case, which involves non trivial cancellations 
between different Feynman diagrams, for the process qq — > scalar unparticles. 



1 Introduction 



Recently Georgi [1, 2] has suggested a new approach to analyzing the phenomenology of 
conformal sectors coupled to the standard model by non-renormalizable interactions. He 
refers to the propagating degrees of freedom of such conformal sectors as "unparticles" since 
their phase space resembles the phase space of a non-integer number of particles. Georgi [1] 
also speculated that it would be interesting to study the case where the unparticles had 
standard model gauge quantum numbers. Since the presence of new massless modes with 
standard model gauge couplings would drastically modify low energy phenomenology, we 
will assume in this paper that such charged unparticles have an infrared cutoff, so that the 
effects of the unparticles will be mainly restricted to high-energies such as will be probed by 
the Large Hadron Collider (LHC). 

The most spectacular signals for unparticles at the LHC would arise in the case where 
the new unparticles carry a color charge. For example we could imagine a theory with a new 
type of quark that carries color but also couples strongly to an approximate conformal field 
theory (CFT) sector. If sufficiently stable, these "unquarks" would appear in the detector 
as "unjets" which can be very different from the usual jets of QCD. Since, as mentioned 
above, we do not think there can be new light degrees of freedom that couple to gluons, the 
CFT must have some dynamics that cuts off the infrared at least below the few hundred 
GeV range. This is similar to the the case of Randall- Sundrum models [3] with gluons in the 
bulk: an infrared cutoff is required to have a phcnomcnologically acceptable model, and the 
usual infrared cutoff (a TeV brane) violently breaks the approximate conformal symmetry 
of the model. We will similarly impose a hard infrared cutoff on the colored unparticles, but 
it is still possible that color neutral unparticles can be massless or at least much lighter than 
the colored unparticles. 

Essentially we are thinking of an effective theory for a massive particle that couples to a 
CFT so that it has a large anomalous dimension. One of the intriguing features of Gcorgi's 
proposal is that it suggests new ways to impose an infrared cutoff on a CFT (or a Randall- 
Sundrum-like) theory. Further it may provide a weak coupling description of a CFT without 
relying on a large N limit. 

Here we will focus on an analysis of colored scalar unparticle production from a quark- 
antiquark initial state. Gluon fusion will probably dominate at the LHC, but that calculation 
is much more involved and we will be content here with demonstrating the general methods 
that go into such calculations. We will nevertheless provide a simple answer for the gluon 
fusion case as well. After introducing our infrared cutoff, we examine colored unparticle 
propagators and how to derive their gluon coupling Feynman vertices (form-factors). Fi- 
nally we describe the parton level production cross-section calculation and finish off with 
conclusions and a discussion of future directions. 

2 An Infrared Cutoff 

The phase space factor for a scalar unparticle final state with scaling dimension d is [1] 
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d^p,d) = A,9{p')9{p') {p'Y-\ (2.1) 

where 

_ 167r5/^ r(rf+l/2) 

(27r)2<i r(t/- l)r(2c/) ^ • ' 

amounts to a useful normalization convention [1]. 

Without knowing the details of how the scale invariance of the CFT sector is broken 
in the infrared, we will parameterize our ignorance with an infrared cutoff scale m. For a 
Lorentz scalar operator with scaling dimensior0 1 < d < 2 we take the simple ansatz [4] : 

A{p,m,d) = [ d'^xe'P^0\TO{x)O\0)\0) 



This form has been chosen because in the limit m — > it reduces to Georgi's unparticle 
propagator, and furthermore in the ci — > 1 limit it reduces to a free particle propagator 
with mass m, which we know to be the case for gauge invariant operators from unitarity 
arguments [5]. 

From the conjectured correspondence between five dimensional anti-de Sitter (AdSs) 
space theories and CFTs, it also seems reasonable to suggest that the infrared cutoff should 
not correspond to something like confinement since this would just result in towers of CFT 
"mesons" , "baryons" , and "glueballs" , which are the Kaluza-Klein modes of the AdSs theory. 
Perhaps the CFT we are considering here has "chiral symmetry breaking without confine- 
ment". There are well understood CFTs where we can engineer such behavior. Consider a 
Banks-Zaks fixed point theory [6] with a large N gauge group. Giving a mass to one of the 
fermions will give an order 1/N correction to the fixed point, but the massive fermion will 
still be able to emit massless gauge bosons and will acquire a branch-cut in its propagator 
starting at the fermion mass. By looking at supersymmetric examples we can extend this 
analysis to strong coupling. Consider Seiberg's dual of SUSY QCD [7] in the conformal 
window at large A^. The meson operator has a scaling dimension between 1 and 2. We can 
assume for simplicity that the gauge coupling is tuned to its fixed point value in the UV. 
Adding a mass term for one of the flavors will again give an order 1/A^ correction to the fixed 
point and result in a propagator with a branch cut starting at the mass. Such branch cuts 
are ubiquitous in field theory. They occur in any theory where a particle can emit massless 
states such gauge bosons. This happens for the electron propagator in QED [8] and the 
gluon propagator in QCD [9], for example. 

Performing the integration in (12.31) we find 

A{p,m,d) = ,r^— ! (2.4) 



2 sin dir {jrfi — p^ 



^For scaling dimensions d > 2 we need to provide an ultraviolet cutoff as well. 
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In the limit d 1 we have 

A(p,m,l) = ^y— -, (2.5) 

the free particle propagator, as required by unitarity [5]. 

The discontinuity across the cut gives the modified phase space to be: 

d<l>{p, m, d) = Ade e _ m^) [p^ - m^) "'"^ . (2.6) 

As d — > 1 from above, d(^{p,m,d) approaches 1-particle phase space with the appropriate 
mass-shell constraint 

d^{p, m, 1) = 27r (/) 5{p^ - m^) . (2.7) 



3 Gauge Interactions 

The propagator in Eq. (12.41) can be derived from the momentum space effective action 



S 



2 sin c/tt 



d^p 



(3.1) 



If we want to weakly gauge a global symmetry G of this theory we can first Fourier Transform 
the action (13. ip to arrive at a non-local theory in position space 



S 



d'^x d'^y (j)\y)F{x — y)4){x) , 



(3.2) 



and then ensure G gauge invariance by using Mandelstam's method [10] of introducing a 
path-ordered exponential of the gauge field, i.e. a Wilson fine. 



W{x, y) = Pexp 



(3.3) 



between the two unparticle fields evaluated at x and y. This method has been previously 
applied to non-local toy-model field theories [11,12], for example the non-local chiral-quark 
model [12]. 

Applying this method to the QCD interactions of the unparticles allows us to calculate 
the Feynman vertex for a gluon coupled to two unparticles. The result using Eq. (13. ip is 



tgT^''{p,q) 



6A''''{q)6(f)^p + q)6(f){p) 



igT 



2sinc?7r 2p^ + q^' 



^^2_^p^^yY-'-{m^-p^y . (3.4) 



Ad 2p- q + q^ 

As a check we note that this vertex satisfies the Ward-Takahashi identity [13] 

iq^V^" = A-\p + q, m, d)T'' - T^A~\p, m, d) . (3.5) 
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The path-ordered exponential includes arbitrarily high powers of the gauge field, so there 
are vertices with arbitrary numbers of gauge bosons. The two gauge boson vertex is 



{2p + q2n2p + 2q2 + qi 



1^ 



qf + 2{p + ga) ■ qi 
+ qiY{2p + 2qi + : 
ql + 2{p + qi) ■ q2 



2 



where 



(3.6) 



Ad q^ + 2p-q 

4 Unparticle Production 

In this section we want to compute the cross section cr{qq unparticles) at leading order in 
as- Thus we need to evaluate the imaginary part of qq forward scattering via a gluon in the 
s-channel with a one loop correction to the gluon line from unparticles, as shown in Fig. [H 
A slick way to get the result is to calculate the effective action by taking the logarithm 
of the partition function Z. From writing fl3.2l) as a power of a gauge covariant derivative 
D, we simply find from the path integral that: 

InZ = -^\nDet{D^ + my-^ 
= -^Trln(D2 + ^2)2-d 

= --(2-ci)Trln(D2 + ^2^)_ 
2 

Thus the sum of one unparticle loop graphs with a fixed number of gauge boson legs is just 
2 — d times the usual one particle loop answer. 

We will verify this result by explicitly calculating the cross section in the qq channel. One 
might think of applying the usual Cutkosky rules [14]: cutting the unparticle propagators and 
replacing them with the unparticle phase space factor (12.61) . This however leads to the wrong 
answer, in particular a divergent cross section. The reason is that the vertices in Eqs. (13.41) 
and (13. 6p also contain non-trivial functions of the momenta, and they do contribute to the 
imaginary part of the qq scattering amplitude. Therefore, we must consider the amplitude 
as a whole and compute the imaginary part directly. 

As a starting point, we can simplify the calculation by setting the infrared cut-off to 
zero. The imaginary part can be extracted from the vacuum polarization amplitude: for 
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Figure 1: Forward scattering of qq through a gluon with a scalar unparticle vacuum polar- 
ization loop. The imaginary parts of those diagrams contribute to the unparticle production 
cross-section. 



incoming momentum g, gauge invariance requires that such an amplitude is proportional to 
g^(ln(— g^/A^) + finite), and the imaginary part can be extracted from the In term 



Im [g2in(-g7A2)] = . 



(4.2) 



Since dimensional analysis requires that the ln(— g^) comes with a logarithm of the regulator 
scale squared, A^, we can easily find the coefficient of ln(— g^) by extracting the regulator 
dependence. 

The first diagram in Fig. [1] is proportional to: 



d^p q^ — Ap ■ q + Ap^ 

(27r)4 (g2_2p.g)2 



— {p — qY — 
— p2 — ie 



2-d 



-p — le 



— [p — qY — ie 



2~d 



(4.3) 



Note that the normalization Ad cancels between propagators and vertices in the loop as it 
should. Looking at the large p region of the loop integration we can Taylor expand for small 
g: 



d'^p q^ — Ap ■ q + Ap'^ 
(27r)4 (g2-2p-g)' 

where 



C2- 



(g^ — 2p ■ qf (g2 — 2p ■ qf (g^ — 2p ■ qf 



p" 



+ C3- 



+ C4- 



+ 



,(4.4) 



C2 
C3 
C4 



id 



{d-2)\ 

-{d-2f (d^ -Ad +15) /12 . 
Picking out the terms that generate logarithms we are left with 



(4.5) 



d'^p , . X g^ / N8(io-g)' 



(4,6) 



Thus, the first diagram's contribution to the massless scalar unparticle production cross- 
section is suppressed relative to the massless scalar particle production by 



"A 



d{2-df{A-d) 



diag. 1 



(4.7) 



The second diagram does not give any contribution in the particle case {d = 1), where 
the only contribution to the imaginary part comes from cuts of propagators. However, this 
is not true in the unparticle case. The amplitude is proportional to the integral: 



A, 



(2^ 



2-d 



2 — m? + ie 



— Ap ■ q + Ap^ ^ q'^ + Ap ■ q + Ap^ 



(g2 _ 2p ■ g) 



(g2 + 2p ■ g) 



+ 



+ 



q"^ — Ap ■ q + Ap^ 
(g2 - 2p ■ qf 

q^ + Ap- q + Ap^ 



[m? — {p — qY — ie 
[m? — p"^ — ief'~'^ 



,2-d 



{w? — {p + qY — ie 



2-c 



9 9 • \2— d 

— p^ ~ le) 



(g2 + 2p ■ qY 
Applying the same technique as before, we find: 



(4.8) 



As expected, this contribution to the cross section vanishes both in the d 1 and d 2 
limits. Combining the two results in Eqs. (14.71) and (14. 9 p together, we obtain the simple 
result 



{d-l)id-2){d' -5d+3) 



diasr. 2 



(4.9) 



- = {2-d), (4.10) 
0-1 

in agreement with Eq. (14.11) . 

In the Appendix we show how to obtain the same results via a direct calculation of the 
discontinuity of the integrals, providing a check of the results in Eqs ( 14.71) and (14. 9p . Such 
a calculation allows us to take into account the presence of the IR cutoff, which cannot be 
inferred from the small q expansion above. Our simple ansatz for the breaking of confor- 
mal invariance, leads to a very simple result: the d and m dependence factorize, and the 
dependence on the IR cutoff is given by the usual kinematical factor. This result is achieved 
via non trivial cancellations between the two diagrams in Fig. [T], as shown in the Appendix. 
Therefore the unparticle production cross section is: 




ad{m) = (2 - rf) Wl - ^ ai(m = 0) = {2 - d) a,{m) . (4.11) 



5 Conclusions 

In this paper we have considered scalar unparticles that carry standard model gauge quan- 
tum numbers. We have seen that the production cross section of such unparticles via gauge 
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interactions is, at first order in Ug, simply suppressed by a factor 2 — d compared to the pair 
production of the corresponding scalar particles. We explicitly calculated the cross section 
in the qq annihilation channel, where this simple result arises through a non trivial cancella- 
tion between two diagrams, one of which is does not contribute in the particle limit d ^ 1. 
We included a simple parameterization of an IR cutoff, which seems phenomenologically 
necessary for unparticles to be charged under standard model gauge groups. Our parame- 
terization yields a simple result: the IR cutoff dependence factorizes and the cross section 
is multiplied by a simple suppression factor. These results may be generalized to the gluon 
fusion channel, which will dominate for unparticle pair production at the LHC. It would be 
quite interesting to study the detailed phenomenology of colored unparticles at the LHC. 

In addition to having a suppressed cross-section, the final states arising from these colored 
unparticles would look very different from ordinary jets. If the unparticles are stable enough 
to travel a significant distance through the detector one would have to understand the details 
of unparticle QCD hadronization ("unjet" formation). Some of the difficulties involved are 
analogous to those of the long-lived gluino scenario [15]. A further comphcation is that the 
unquarks can still radiate into hghter CFT degrees of freedom that are not colored (the 
gauge singlet unparticles usually considered in the literature). Thus there is an additional 
loss of energy that is invisible to the detector, which results in missing energy along the 
two unjet directions. In a more general model than we have considered here, with mixing 
between colored unparticles and other colored particles, one could imagine a scenario where 
a single colored unparticle is produced so that the missing energy is entirely along a single 
unjet direction. 

An even more interesting scenario to explore is one where the unparticles also have 
electroweak quantum numbers. Unquarks could then decay into ordinary quarks and they 
would have a passing resemblance to a fourth family [16]. (This is the simplest way to allow 
unquarks to decay and avoid searches bounding stable exotics.) We could even imagine 
that the unquarks play a role in electroweak symmetry breaking which would give us a new 
approach to the hierarchy problem. If the scaling dimension of the analog of the Higgs mass 
operator is larger than two, such a scenario would have some similarities with gaugephobic 
Higgs models [17]. 
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A Appendix: Calculation 



In this appendix we discuss an alternative way to calculate the integrals in Eqs (14.31) and 
(14. Sp . that correctly takes into account the effect of the IR cutoff of the conformal sector. 
Let us focus on Eq. (14.31) : 



d^p q"^ — Ap ■ q + Ap^ 
(27r)4 {q^-2p-qf 



1 1 
m ~ p — le 

rn? — {p — qY — ie 



2-d 



— {p — qY — ie 
m? — p"^ — ie 



2~d 



(A.i: 



We are interested in the imaginary part of ( lA.ll) . so we can neglect the factor 2 inside the 
square brackets. For simplicity, we will illustrate the calculation on the second term in the 
integral: 



d'^p q'^ ~ Ap ■ q + Ap^ 



1 1 
m — p — le 



(27r)4 (g2_2p.g)2 - - q)2 _ 

We use Feynman parameters to rewrite this integral as 



2-d 



(A.2) 



- sinrfvr 
Im :i[q) = Im 



TT 



d^p q'^ — Ap ■ q + 4p 
(27r)4 (g2_2j9-g)' 



2 /"I 



dx 



'^(1 — xY ^ (p^ — m^) 
p2 + x(g2 — 2p ■ q) — w? + ie 



(.A.3) 



In this way the rf— dependence appears only in the powers of x: momenta all have integer 
powers. Going to the rest frame, where p = {E, p),q = {q, ), we notice that the integrand 
has a pole at E = qx + ^Jp^ — — x)x. Taking the residue at this pole we get 



Im H(g) = — ^ ^ ^ Im / dp j dx ^ 



1 — X 



2-d 



p^ {Am? + q {A^Jra^ + p^ — q^{l — x)x{2x — 1) + g(l — 8(1 — x)x) 



(A.4) 



q^J m? + p^ — g2(l — x)x {q{2x — 1) + 2^Jm? + p^ — q'^{l — x)x^ 
Next notice that the term (g(2x — 1) + 2^vrfi + — cp'iX — appearing at the denomi 



nator of Eq. (1A.4I) vanishes at p = \/ q^ — Am? /2 only if a; < 1/2. Thus this pole will give a 
contribution to the imaginary part of the integral. We evaluate it taking the residue 



Impoic S(g) 



sin (ivr / dx 



327r2 

J2 



1 — X 



^ sin dn p-" B { -; 3 - rf, -1 + 
327r'' 



(A.5) 



where f3 = a/1 — Am^ /q^ and B is the incomplete beta function, defined as 
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B{z;a,b) = / du {uf-\l-u) 
Jo 



b-l 



(A.6) 



There is another contribution to the imaginary part of our original integral, which comes 



directly from the region where the argument of the square roots appearing in Eq. (IA.4p is 
negative. This can be evaluated as 



Im^ E{q) 



ax 



^y—m'^+q^x(l—x) 



1 ~ X 



2~d 



dp 



{Am? + - qi)(l- 2x) 



(4 [w? + p^) — q^) ^J —vrfi — + — x)x 



(A.7) 



The integral can be performed analytically and the result can be expressed in terms of 
incomplete Beta functions: 



Im^ 



+ 



sm avr 

327r2 

/33 + 2/32-1 



5 ( ^; 3 - rf, d - 1 ) + /^^ Bi}-^ 



■,3 — d,d — 1) + 



fiP, 3 - rf, d - 1) + (3 - 2p^) f{p, 4-d,d-l) + 



6 5 - rf, - 1) + 4 /(/?, 6 - rf, - 1^ 



where we defined the auxiliary function 

fiP,a,b) : 



B{^—^;a,b] . 



(A.9) 



It is interesting to note that the contribution of the last term in Eq. lA.ll can be obtained by 
the previous result substituting d ^ 4 — d. 

The contribution to the cross section of the 2 diagrams in Fig. [1] can be written as: 



(Td(m) 



diag. 1 



sin dn 



71 



+ 2/3^-1 



(/(A 3-d,d-l)- /(/?, d-l,3-d)^ + 



+ P^B{- — ^3-d,d-l)~B{- — ^■,d-l,3-d)] + 



+ (3-2/?2)(/(/3,4-rf,d-l)-/(Af/,3-rf)) + 
- 6(^/(/3,5-rf,rf-l)-/(/3,d+l,3-d)) + 

+ A(^f{(3,Q-d,d-l)~ f{(3,d + 2,3-d) 



(A.IO) 
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diag. 2 ^1 



(A.ll) 

diag. 1 



where the result from the second diagram matches non trivially the same series of Beta 
functions as the first diagram, leading to the simple result presented in Section HI 

We can now check that in the case m — we reproduce the results in Eqs (14. 7p and 
(14. 9p . using the following properties of the Beta functions: 

B(0;a,b) = 0; 5(1; a, 6) = ^^^^^ , if a,b>0. (A.12) 

r(a + b) 
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